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Abstract 
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A Fourier inversion problem lies at the heart of determining spatio-temporal 
characteristica of the final stage of a heavy ion collision: From the measured 
two-particle momentum correlations C(pi,p2) of identical particles, pions say, 
a Hanbury-Brown /Twiss (HBT) interferometric analysis aims at extracting as 
much information as possible about the Wigner phase space density S(x,p) of 
pion emitting sources in the collision region, Jl], |2| 
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where q = Pi — P2, K = \{jPi+P2)- Here, we discuss how the analysis of (|l|) 
allows to separate the effects of temperature and transverse flow which cannot 
be disentangled completely on the basis of single-particle spectra. 
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1 Paper presented at the XXXIIth Rencontres de Moriond, "QCD and High Energy Hadronic Interac- 
tions", Les Arcs, France, March 22-29, 1997. 



Physics is rich in Fourier inversion problems which are not analytically invertible. A 
well-known example is the determination of the charge density distribution p(r) of atomic 
nuclei from the structure function F(q) measured in elastic electron scattering, 



F 2 (q) = |i/rfrp(r)e^| 2 . (2) 

In most cases, this structure function was measured only in a small window of the momen- 
tum transfer q which makes an analytical unfolding of p(r) without additional assumptions 
impossible. The standard way out is to take recourse to model distributions p(r) like e.g. 
the Fermi distribution p(r) = po/ (l + exp( L j^)^J whose parameters have a simple physical 
interpretation. In a first step, these model parameters are determined in a fit of F(q), till 
a refined analysis reveals the need for more detailed models which include e.g. dips in the 
charge distribution for small r, cf. ||. 

In our analysis of the measured two-particle correlations ([I]), we adopt essentially the 
same strategy, for different reasons though: the detected pions which determine C(q, K) 
are on mass-shell and this results in a constraint for the relative pair momentum q, q-K = 0. 
Hence, the four-vector q has only three independent components and even for an exactly 
known correlator, the Fourier transform (|XJ) is not uniquely invertible. This necessitates a 
model-dependent approach. 

The analysis presented here is based on a model pion emission function including reso- 
nance decay channels R: 

S„(x, p) = S^ T (x, p) + J2 S R->A X , P) ■ ( 3 ) 

R 

The contributions are obtained from the direct resonance emission function Sr(X, P) 
by propagating the resonances of widths T, produced at (X^,P^), along a classical path 
x^ = X M + j^-t according to an exponential decay law j[| 

Sr^(x, P ) = j^j d 4 X J drTe- rT S^(x - (x + ^r))Sf(X,P) , (4) 

where J* R is the integral over the available resonance phase space for isotropic decays. 
Our model assumes local thermalization at freeze-out and produces hadronic resonances 
by thermal excitation. For particle species % with spin degeneracy 2Jj + 1, the emission 
function reads || 

Sf^, P) = P-<x) exp exp (- £ - ^ - . (5) 

The Boltzmann factor exp[— (P-u(x) — fii)/T] implements both the assumption of ther- 
malization, with temperature T and chemical potential pj, and collective expansion with 
hydrodynamic flow 4- velocity u^(x). The geometrical extension of the collision region is 
determined by Gaussian widths R and A77 in the transverse radius r and the space-time 
rapidity 77 = |ln [(t + z)/(t — z)], as well as a Gaussian average around a mean freeze- 
out proper time r with dispersion At. Freeze-out occurs along proper time hyperbolas 
P-n(x) = Mj_ cosh(y — if) where Y and Mj_ are the rapidity and transverse mass asso- 
ciated with P. In the longitudinal direction we assume scaling expansion, v\ = z/t or 
i]i = |ln[(l + ^)/(l — v{)] — T), for the transverse expansion a linear rapidity profile: 




We include all pion decay channels of p, A, K*, £*, uj, p, p', K s , S and A with branching 
ratios larger than 5 percent. 

For a given model emission function, both the two-particle spectrum (fj) and the one- 
particle spectrum E p dN/d 3 p = J d 4 x S(x,p) is determined. On the LHS of Fig. [I], we plot 
the rapidity (y) integrated pion transverse mass spectrum 
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for two sets of source parameters without (pj = 0) and with (pj = 0.3) transverse flow. 
Note that the geometric parameters R, Ap, r of the source enter only in the normaliza- 





Figure 1: Left: One pion transverse momentum spectra according to (0) for different 
values of the transverse flow pf. Right: The cr-confidence levels for a fit of our model to 
the NA49 /i~-spectrum at midrapidity. 

tion of this spectrum. The slope is essentially determined by the temperature T and the 
transverse flow pf and contains no information on the source geometry. The direct pions 

reflect essentially an effective "blueshifted" temperature T e s = Ty ^^ ". If resonance 

contributions to (0) are included, these statements change in many details, most notably, 
resonance contributions lead to a low pr-enhancement, and a more pronounced flattening 
of the spectra for larger transverse flow. The main message however remains unaltered: 
the unnormalized one-particle spectra (0) contain no geometrical information, their slope 
determines only a combination of temperature and transverse flow. This is clearly seen on 
the RHS of Fig. p] where we show the confidence levels of a fit to the NA49 /^"-spectrum 

Two-particle correlations are often analyzed in terms of the Gaussian ansatz 

C(q, K) = 1 + A e~^~ R ° q °~ R ^ q t~^ q ° qi . (8) 

They allow to disentangle temperature and flow effects. To see this, one can consider e.g. 
the approximate expression 
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Here, the m_i_-dependence of (||) is governed by a prefactor ifj/T. This indicates that the 
side radius R S (K±) may allow to distinguish scenarios with a relative large temperature and 
a small transverse flow size rjf from those with a large rjf and small T which according to 
the RHS of Fig. [I] account for realistic transverse momentum slopes equally well. We hasten 
to remark however that (||) is obtained in a crude saddle point approximation neglecting 
all resonance decay contributions - it is known to be quantitatively unreliable. This was 
one of our main motivations for a quantitatively accurate, numerical calculation || of the 
two-particle correlation C(q, K) and its Gaussian widths. 

Most notably, resonance decay contributions can lead to deviations of the correlator 
from a Gaussian shape. This is a consequence of the exponential resonance decay law 
and can lead to significant quantitative ambiguities in analyses based on (|8]), see e.g. the 
different Gaussian fits depicted in Fig. |2|. We hence suggest || to characterize the correlator 
by a set of so-called g-variances which do not depend on any assumption about the shape. 
For a unidirectional analysis of the correlations C( qi , K) = C( qi , g,,yj=0, K) along the three 
Cartesian axes, these g-moments read 
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Note that for a correlator of Gaussian shape, these expressions coincide with those of a fit 
to (|8|). Still, they are well-defined for arbitrary shapes of the correlator. Deviations from 
a Gaussian shape can be discussed in terms of higher g-moments, most easily in terms of 
the "kurtosis" 

^SW" 1 ' (13) 

Fig. ^| shows the numerically obtained results for the side radius if resonance decay contri- 
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Figure 2: HBT radius parameters R s and the kurtosis A s calculated from the g- moments 
and compared to a Gaussian fit to (§). The two plots depict scenarios without (rjf = 0) 
and with (rjf = 0.3) transverse flow. 

butions are added: R s develops a -fd-slope even for scenarios without transverse flow. The 
reason is that the resonances propagate outside of the thermally equilibrated region before 
decaying. This "lifetime effect" leads to exponential tails in and tends to increase the 

Gaussian widths Ri. Since the relative abundance of resonances is larger for small K±, this 



tail is more pronounced in the region of small K±, thereby leading to a Xj_-slope of R s . In 
contrast, for a finite transverse flow 77/ = 0.3 the i^j_-dependence of R s is due to the flow, 
and resonance decay contributions do not change the slope of R s in our model. This can 
be traced back to a shrinking transverse size of the direct resonance emission function Sf{ r 
which for the case depicted in Fig. ||, counterbalances the lifetime effect almost exactly |J . 

Here, the main message is that resonance decay contributions can lead to a i^-slope 
of R s , even in scenarios without transverse flow, i.e., they are faking a signal previously 
attributed to transverse collective dynamics. However, the physical origin of the _?G_-slope 
of R s is different for the two scenarios in Fig. |2| and it is well understood || that this 
physical difference manifests itself in a much more prominent "non-gaussicity" for the case 
7]f = 0. The kurtosis ( |T3"D allows for a quantitative measure of deviations from a Gaussian 
shape and hence it provides at least for the class of models discussed here, the cleanest 
distinction between scenarios with and without transverse flow. We hence expect that in a 
detailed analysis of two-particle correlation functions, based on Eqs. (|i~0"l) - (|i~3| ), the effects 
of temperature and transverse flow can be disentangled most clearly. 

In this short presentation, we have focussed entirely on one aspect, the distinction 
between thermal excitation and transverse collective dynamics in the collision region. We 
conclude by pointing out that various other geometrical and dynamical characteristics of 
the source are accessible via the analysis of two-particle correlations. Especially, we mention 
information on the longitudinal expansion, emission time and emission duration as well as 
the longitudinal and transverse size of the boson emitting region. In the context of the 
model (|5|), resonance decay contributions to these characteristica are quantitatively under 
control and we are currently working on the determination of the model parameters of @ 
from the data [|?J taken by the NA49 Collaboration. 

This report is based on joint work with Ulrich Heinz. It was supported by BMBF, 
DFG, GSI and DPNC. 
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